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ABSTRACT 

We estimate the fraction of mass that is composed of compact objects in 
gravitational lens galaxies. This study is based on microlensing measurements 
(obtained from the literature) of a sample of 29 quasar image pairs seen through 
20 lens galaxies. We determine the baseline for no microlensing magnification 
between two images from the ratios of emission line fluxes. Relative to this 
baseline, the ratio between the continua of the two images gives the difference in 
microlensing magnification. The histogram of observed microlensing events peaks 
close to no magnification and is concentrated below 0.6 magnitudes, although two 
events of high magnification, Am ~ 1.5, are also present. We study the likelihood 
of the microlensing measurements using frequency distributions obtained from 
simulated microlensing magnification maps for different values of the fraction of 
mass in compact objects, a. The concentration of microlensing measurements 
close to Am ~ can be explained only by simulations corresponding to very 
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low values of a (10% or less). A maximum likelihood test yields a = 0.05 
(90% confidence interval) for a quasar continuum source of intrinsic size r So ~ 
2.6 ■ 10 15 cm. This estimate is valid in the 0.1 — 1OM range of microlens masses. 
We study the dependence of the estimate of a with r so , and find that a < 0.1 
for r S[) < 1.3 • 10 16 cm. High values of a are possible only for source sizes much 
larger than commonly expected (r So >> 2.6 • 10 16 cm). Regarding the current 
controversy about Milky Way/LMC and M31 microlensing studies, our work 
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supports the hypothesis of a very low content in MACHOS (Massive Compact 
Halo Objects). In fact, according to our study, quasar microlensing probably 
arises from the normal star populations of lens galaxies and there is no statistical 
evidence for MACHOS in the dark halos. 

Subject headings: gravitational lensing, dark matter, galaxies: halos 



1. Introduction 



The composition of matter in the halos of galaxies is a central problem in astrophysics. 
Duri ng the last 10 y ears, several observational projects have used gravitational microlens- 



ing (IPaczvn ski 



Alcocketal 



19861 ) to probe the properties of the halos of the Milky Way (MACHO, 



2000 ; EROS, 



Tisserand et al.ll2007h andM31 (POINT- AGAPE. ICalchi Novati et al 



20051 : MEGA, Ide Jong et al.ll2006l ). These experiments are based on the detection of mag- 
nification in the light-curve of a source induced by an isolated point-like (or binary) object 
passing near the observer's line of sight. From the successful detection of a number of mi- 
crolensing events these collaborations have estimated the fraction of the halo mass that is 
composed of lensing objects, a. However, the reported r esults disagree. For the Milky Way's 
halo the measurements of the MACHO collabor ation (lAlcock et al.l 120001) correspond to a 
halo fraction of 0.08 < a < 0.50 while EROS JTisserand et al.l |2007h o btains a < 0.08. 



On t he other hand, re-analysis of publicly available MACHO light-curves (IBelokurov et al. 



20041) leads to results simila r to those reported by E ROS (however, see also the counter-report 
bv lGriest fc Thomasl I2OO5} ) . For M31 the AGAPE JCalchi Novati et al.l l2005h collab oration 
finds a halo fraction in the range 0.2 < a < 0.9, while MEGA ( Ide Jong et al.l 120061 ) finds a 
limit of a < 0.3. 

The method applied to the Milky Way and M31 can be extended to the extragalac- 
tic domain by observing the microlensing induced by compact ob jects in the lens galaxy 
halo in images of mu ltiply imaged quas ars (quasar microlensing; IChang fc Refsdall Il979l . 
see also the review by IWambsganssI 120061 ). Interpreting the light-curves of QSO 2237+0305, 
Webster et al.l (Il99ll ) suggest that the monitoring of microlensing variability can provide 
a measure of the op tical depth in compact objects and in the smooth mass distribution. 
Lewis fc Irwinl (119961 ) proposed a statistical approach to the determination of the mass den- 
sity in compact objects based on the comparison between the observed and simulated mag- 
nification probability distributions. Microlensing can also be measured from a single-epoch 
snapsho t of the anomalo us flux ratios induced by this effect bet w een the images of a lensed 
quasa r (IWitt et al.lll995l : see also lSchechter fc Wambsganssll2002l ). ISchechter fc Wambsganss 



(120041 ) explore the practical application of this idea by using a sample of eleven systems with 



- 3- 



measured flux anomalies. Other quasar microlensing studies of interest for the present study 
are aimed at the determination of accretion disk sizes (e.g. , the studies based in relatively 



large samples by iPoindexter et al.l 120071 . iMorgan et al.l 120071 and references therein) . 



In practi ce, the study of extragalactic microlensing meets significant obstacles, in par- 



ticular (e.g., iKochanek! 12004 ) larger time-scales for microlensing variability and lack of a 
baseline for no magnification needed to detect and to quantify mi crolensing (see, howe ver, 
the time variability based studies of several individual systems in IMorgan et al.l 120081 and 
references therein). In addition, microlensing by an isolated object is not a valid approxima- 
tion. Mic rolensing at high optica l depth should be modeled (e.g., by simulating magnification 
maps; see 



Schneider et al.lll992f ). 



We avoid these obstacles by setting the baseline of no microlensin g magnification using 



the n arrow emission lines (NELs) in the spectra of lensed quasar images (jSchechter fc Wambsganss 



2004 follow a similar approach but using theoretical models to define the baseline). It is gen- 



erally expected that the regions where NELs originate are very large (compared with the con- 



Kaspi et al. 


2000; 


Aba j as et al. 


— , „ 
2002) 



If we define the baseline from emision lines measured in the same wavelength regions as 
the continua affected by microlensing, we can also remove the extinction and isolate the 
microlensing effects. 

"Intrinsic" flux ratios between the images in the absence of microlensing can be deter- 
mined from the observation of the mid-infrared and radio-emitting regions of quasars that 
should also be large enough to average out the effects of microlensing (see IKochanek! 12004 
and references therein). However, the extinction at mid- infrared and radio wavelengths is 
lower than the extinction at the wavelengths in which microlensing is usually detected and 
measured (optical, near-infrared, and X-ray). Consequently, the difference between the mid- 
infrared (radio) and the optical (X-ray or near-infrared) continuum fluxes will include not 
only the effects of microlensing but also the effects of extinction. In addition, note that the 
availability of data at optical wavelengths is considerably greater than at other wavelengths. 

Thus we will use the NEL and continuum flux ratios among the different images of a 
lensed QSO to estimate the difference of microlensing magnification between the images at 
a given epoch with certain restrictions that we detail in the following paragraphs. 

The flux (in magnitudes) of an emission line observed at wavelength A of image i of a 
multiply imaged quasar is equal to the flux of the source, mj)" 1 > magnified by the lens 

galaxy (with a $j magnification factor; /ij = —2.5 log $j) and corrected by the extinction of 
this image caused by the lens galaxy, Ai 



(see, e.g.. lMunoz et al.ll2004l ) 
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m? n (\) = m% 



tin 



l + Zs 



+ fa + Ai 



l + zi 



where z s and zi are the redshifts of the source and the lens, respectively. 

In the case of the continuum emission, we must also take into account the intrinsic 
variability of the source combined with the delay in the arrival of the signal, which is different 
for each image, Atj, and the microlensing magnification, which depends on wavelength and 



time (with a 



1 + 2; ' 



magnification factor; A/Xj = — 2. 5 log < 



mr(A,t) = m- 



A 



l + z. 



t- Ati) + fa + Ai 



l + zi 



i + V 



t). 



(2) 



Thus, the difference between continuum and line fluxes cancels the terms corresponding 
to intrinsic magnification and extinction (/jj + Ai): 



mf» (A, t) - m l r{\) = mS m (j^-, t - At^j - (r^A + (j^> *) ■ ( 3 ) 

If we consider a pair of images, 1 and 2, the continuum ratio relative to the zero point 
defined by the emission line ratio can be written (in magnitudes) as, 

Am(A, t) = {m x - m 2 ) con - {m x - m 2 ) lm = A/ii (r^> *) ~ A ^ (i^, *) (4) 

+Amr (^^tx-Ata). 

We have referred the equations for the magnification of both images to an arbitrary 
time, t (note that microlensing-induced variability between a pair of images is uncorrelated) . 

The first term of Equation (j3J) is the relative microlensing magnification between images 
1 and 2. The significance of the second term, Am™", which represents the source variability, 
can be estimated by comparing the intrinsic quasar variability on time-scales typical of the 
time delay between images in gravitational lens systems with the expected distribution of 
microlensing magnifications. As we shall discuss in Section HJ the intrinsic source variability 
is not significant for our computations. 

In summary, with the proposed method similar information as in the Milky Way MA- 
CHO experiments is obtained but with a single-epoch measurement. The objective of this 
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study is to apply this method to published data of quasar microlensing. In $2] we collect the 
data from the literature and fit models to the systems of multiply imaged quasars to derive 
suitable values of the projected matter density at the image locations. Using these values, 
probabilistic models for microlensing magnifications are derived in Section 13.11 for a range 
of fractions of mass in compact objects. Sections 13.21 and 13.31 are devoted to estimate this 
fraction. Finally, in §4] we present and discuss the main conclusions. 



2. Observed Microlensing Magnifications and Macro-Lens Models 

We collected the data, Am (see Eq. @J, examining all the optical spectroscopjQ found 
in the literature (see Tabled]). In most cases the microlensing magnification or the scaling 
of the emission line ratio with respect to the continuum ratio are directly provided by the 
authors or can be estimated from a figure. For SDSS 0806+2006, FBQ 0951+2635, SDSS 
J1001+5027, QSO 1017-207, SDSS J1206+4332, HE 1413+117, and SBS 1520+530 we used 
the electronically available or digitized spec tra of the images to es timate the microlensing 



magnification following the steps described in lMediavilla et al.l (120051 ). In Table[T]we include 



(when available) the flux ratios for each line and its corresponding continuum. Specific details 
of the procedure followed to obtain the data are also given. 

For some of the image pairs (~30% of the sample) there are mid-IR flux ratios available. 
Except for one system, SPS S J1004+4112 (where image C is probably affected by extinction, 



Gomez- Alvarez et al.ll2006l ). they are in very good agreement with the emission-line flux 
ratios (see Table [2]). The average difference between mid-IR and emission line flux ratios 
is only 0.11 magnitudes (0.07 magnitudes if SDSS J1004+4112 is removed). In fact, the 
agreement is unexpectedly good taking into account the possible influence of extinction 
and source variability. In any case, this comparison supports the consistency of the basic 
hypothesis (that the emission line fluxes are not affected by microlensing) and the reliability 
of the data. 

Figure [T] shows the frequency of observed microlensing magnifications, /^(Am). This 
histogram exhibits two significant characteristics: the relatively high number of events with 
low or no microlensing magnification and the concentration (~80%) of the microlensing 
events below |Am| = 0.6. Any model attempting to describe microlensing magnification 



1 There are also sever al X-ray events in the literature that have been explained in terms of microlensing 



(e.g., iPoolev et alJ 120071 and references therein). These events probably arise from a tiny inner region, as 



compared with the optical continuum emitting region, and deserve an analogous but separate study when a 
sufficiently large sample of X-ray microlensing measurements become available. 
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should account for these features. At a lower level of significance, the presence of two events 
of high magnification, Am ~ 1.5, should also be noted. The data presented in Table [1] come 
from many different bibliographic sources with the subsequent lack of information about 
measurement procedures and estimate of uncertainties. However, even with this limitation, 
the low frequency of high magnification microlensing events in the optical seems to be a 
reliable observational result. 

For each of the gravitational lens systems in Table [3] we have used a singular isothermal 
sphere plus external shear model (SIS+7 e ) to estimate values of the total projected matter 
density k and the shear 7 at the locations of the images (see Table [3]). The models have 



been computed with the "lensmodel" code by lKeetonl (1200 ll ) to fit the positions of the images 



(CASTLES, http://www.cfa.harvard.edu/castles/). For double systems we have used the 



emission-line flux ratios between images as an additional constraint. 

In Table H] we show n KltK2 {k\ < k 2 ), the frequency distribution of image pairs that 
occur at combined projected matter densities (ki,k 2 ). The distribution peaks at bin (k\ = 
0.45, k 2 = 0.55). In many of the image pairs in Table [3] the images are roughly located at 
similar distances from the lens galaxy center, 77 ~ r 2 . In an SIS model the convergence for 
each of the lensed images is given by K\ = 1/2(1 + x) and k 2 = 1/2(1 — x), where x is the 
position of the source in units of the Einstein radius. The image configuration n ~ r 2 is 
obtained when x > 0; therefore, the expected values for the convergence are K\ < 0.5 and 
k 2 > 0.5. This is in agreement with Table H] and in fact this simple reasoning could have been 
used to estimate, from a statistical point of view, the peak of the distribution of convergence 
values, n KljK2 , in our sample. 

Observational uncertainties in the flux ratios, differential extinction in the lens galaxy 
and more complicated mass distributions for modeling the lens galaxy could have an im- 
portant impact on the estimates of k and 7. Therefore, the values for k and 7 in Table [3J 
computed for an SIS+7 e model should be taken individually only as compatible values with 
high uncertainties. However, we will assume that the distribution of values in its entirety 
can be considered as statistically representative for the sample of observed image pairs. 
Fortunately, the uncertainty in the macro-lens models does not play a crucial role in the 
conclusions of our study. 



3. Statistical Analysis of the Observed Microlensing Magnifications 



To analyze the microlensing magnification data, we need to consider that each Am mea- 
surement results from the flux ratio of two images seen through different locations at the lens 
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galaxy. The microlensing magnification probability of a given image, //t, >1J K ll 'y 1 (mi), depends 
on the projected matter density in compact objects, the total projected mass density, 
Ki, and the shear, Thus, the probability distribution of the difference in microlensing 
magnification of a pair of images, Am = mi — m 2 , is given by the integral 



(Am) = J / K » llKl , 7l (mi)/ K , 2iK2)72 (mi - Am)dm 1 . (5) 

To simplify the analysis we will suppose that the fraction of matter in compact ob- 
jects, a = k*/k, is the same everywhere. The probability distribution of the difference in 
microlensing magnification of a pair of images can then be written as 



/aKl,Q:K2)Kl,K2,71)72 (Am) J faK\ ,K1 ,71 {j^l ) faK2 ,K2 ,72 (j^l Affljrfffll. (6) 

From this expression we can evaluate the probability of obtaining a microlensing mea- 
surement Am,- from a pair of images, f 1 t 4 < < < .(Am,). Then, to estimate a using 

' c ° ' J CfK^, OK 2 i K li K 2i7li72 

all the available information we maximize the likelihood function corresponding to the N 
measurements collected in Table [Q, 



N 

logL(a) = XI lo s/aKi, Q «*, «-,«*, 7 *, 7 - ( Am <)- ( 7 ) 
1=1 



3.1. Probability Distributions of the Difference in Microlensing 
Magnifications for Image Pairs, / aKl , aK2 , Klift2 , 7l , 72 

We first compute the microlensing magnification probability distributions for one image, 
faK,K,-y{ m )- The first step is to simulate microlensing magnification maps for the different 
values of k and 7 in Table [31 We consider several values for the fraction of mass in compact 
object^: a = 1, 0.5, 0.3, 0.25, 0.2, 0.15, 0.10, 0.05, 0.03, and 0.01. The histogram of each mag- 



2 This sequence of microlensing maps parameterized by a — k*/ K assumes that the overall mass distribu- 
tion (compact o bjects and smooth mass distribution) is close to isothermal. However, in many studies (e.g. 



Dai et al 



2009) the lens galaxy is simulated with a constant mass-to-light (M/L) ratio model representing 



the galaxy stellar content (typically a de Vaucouleurs profile) embedded i n a smooth halo of dark matter 
with no compact objects (a NFW halo, for instance; Navarro et al. 1996h . In this case, the sequence of 
models is parameterized by /m/l, the fraction of mass in the stellar component relative to a constant M/L 
ratio model with no halo (that is, the model with Jm/l = !)• Although the meanings of a and /m/l & re 
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nification map then provides a frequency distribution model of microlensing magnifications. 



We obtain square maps 24 Einstein radii on a side with a spatial resolution of 0.012 
Einstein radii per pixel. To compute the magnificat ion maps we use the inverse polygon 
mapping method described in iMediavilla et al.l (120061 ). An example of the maps is shown in 
Figure [2j The microlensing magnification at a given pixel is then obtained as the ratio of 
the magnification in the pixel to the average magnification. Histograms of these normalized 
maps give the relative frequency of microlensing magnifications for a pixel-size source (see 
some examples in Figure ED . These distributions are in agreement with the results obtained 
with a different method by lLewis fc Irwin! (119951 ). 



To mode l the unresolved quasar source we con sider a Gaussian with r s = 2.6 • 10 15 cm 



(1 light-day) (jShalyapin et al.ll2002l ; lKochanekll2004l ). The convolution of this Gaussian with 
the "pixel" maps gives the magnification maps for the quasar. For a system with redshifts 
z\ ~ 0.5 and z s ~ 2 for the lens and the source respectively, the Einstein radius for a compact 
object of mass M is r] ~ 5.2 • 10 16 a/M/M cm. Thus for M = 1M Q , r] ~ 5.2 • 10 16 cm, and 
the size of a pixel is 6.2 • 10 14 cm. 

Finally, the histograms of the convolved maps give the frequency distributions of mi- 
crolensing magnifications, / aKAi7 (m), which show differences with respect to the results 
obtained for a pixel-size sourc e at t he high magnification wing (the same effect that can 
be observed in iLewis fc Irwinl Il995l ). From the cross-correlation of pairs of these indi- 
vidual probability functions /^^^(mi), and / a K 2 ,K2,7 2 ( m 2) ( see Eq. Mj we obtain the 
probability function of the difference in microlensing magnification between two images, 
/«i,« !lS iA 2 , 11l72 (Am = mi - ra 2 ). In Figures HI and [6] the / aM) „ K2 , K1)K2m , 12 (Am) distri- 
butions corresponding to the 29 image pairs of Table [3] are plotted. 



3.2. Maximum Likelihood Estimate of the Fraction of Mass in Compact 

Objects, a. Confidence Intervals 

In Figure[7Jwe present log L(a) (see Eq. [7]). Using the \ogL(a±na a ) ~ \ogL max —n 2 /2 
criterion we derive a(logL max ) = 0.101q; 03 (90% confidence interval). 

The maximum likelihood method can be affected by errors in the microlensing measure- 
ments, 0"A mi • From eq. [7J we obtain, 



different, the results of both procedures can be compared obtaining from each Jm/l model values of k and 
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JL i df t t * t i A Ami) 
AlogL(a) - \ , fAm .) aAm - (8) 

According to this last expression, microlensing measurement errors do not significantly 
affect the likelihood of flat probability distributions (typical of large values of a). On the 
contrary, the likelihood functions corresponding to low values of a (associated with sharply 
peaked probability distributions) can be strongly modified by the microlensing measurement 
errors. Notice, moreover, that these changes tend to penalize the low a hypothesis. 

To show the impact of (J& mi on the maximum likelihood estimate of a, in Figure [7] we 
also present logL(a) (see Eq. [7j) with error bars, ±AlogL(a), estimated considering that 
each Am; is a normally distributed variable with a^ mi = 0.20 (a realistic estimate). Using 
the logL(a ± na a ) ~ logL max — n 2 /2 criterion and taking into account the error bars of 
logL(a), we derive a(\ogL max ) = 0.05^to!ol (90% confidence interval). 



3.3. Influence of the Continuum Source Size. Influence of the Microlenses 

Mass 

Increasing the size parameter of the Gaussian representing the continuum source, r s , 
affects the previous results by smoothing the magnification patterns and, consequently, the 
probability distributions. To study the dependence of the estimate of a on the source size 
we have computed probability and likelihood functions for several values of this parameter, 
r s = 0.62 ■ 10 15 , 2.6 ■ 10 15 , 8 • 10 15 , and 26 • 10 15 cm. To correct r s from projection effects we 
have taken into account that the intrinsic and projected source areas are related by a cosi 
factor; that is r s ~ \J cos % r so . Assuming that the (disk) sources are randomly oriented in 
space (the probability of finding a disk with inclination, i, proportional to sini) and averaging 
on the inclination, we obtain r so ~ 1.5r s . In Figure [8] we present the likelihood functions 
corresponding to sources of several deprojected size parameters, r So . In Figure M we plot 
the maximum likelihood estimate of a versus^] r So . Error bars correspond to 90% confidence 
intervals. According to this figure, low values of a are expected for continuum source sizes, 
r SQ , of order 1 16 cm or less . Obs erving microlensing variability for nine gravitationally 



lensed quasars iMorgan et al.l (120071 ) measure the accretion disk size. The average value of 
the nine half-light radius determinations is < r\n >= 6 ■ 10 15 cm. For this value we found 



3 Note that for the considered Gaussian intensity profile the radii enclosing 50% and 90% of the source 
energy are related to the Gaussian source size parameter, r So , according to r 1 / 2 = r(50%) = 1.18r so and 
r(90%) = 2.1r so 
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Morgan et al.l (120071 ) report a scaling between the accretion disk 



(see FigureED a = 0.05^;^. 
size an d the black hole mass. In the range of black hole masses considered by iMorgan et al. 
(120071 ) the maximum is Mbh — 2.37 ■ 10 9 M Q , which, using the scaling derived by these 



authors, corresponds to r\j% = 2.4 ■ 10 16 cm. For this size we obtain (see Figure [H]) o. ~ 0.10. 
Values M BH > 10 10 M^ fn/9 > 3.4 • 10 16 cm) should be considered to obtain a > 0.20. On 
the other hand, iPooley et al.l (120071 ) comparing X-ray and optical microlensing in a sample 
of ten lensed quasars inferred r\/2 ~ 1.3 ■ 10 16 cm. For this size we obtain (see Figure [9]) 
a = 0.10to!o6- Thus, according to these recent size estimates based on the observations of 
two relatively large samples of gravitational lenses, high values of a are possible only if the 
continuum source size is substantially larger than expected. 



Owing to the scaling of the Einstein radius with mass, rj oc v M, a change in the mass 
of microlenses can be alternatively seen as a change in the spatial scaling of the magnification 
pattern that leaves invariant the projected mass density, k. Thus, multiplying the mass of 
the microlenses by a factor C (and leaving unaltered the continuum size) is equivalent to 
multiplying the size of the continuum source by a factor 1 / yC (leaving unaltered the masses 
of microlenses). Then the computed models corresponding to sources of sizes r s = 0.62 • 
10 15 , 2.6- 10 15 and 8- 10 15 cm (with 1M microlenses) are equivalent to models corresponding 
to microlens masses of 17, 1, and O.1M (with r s = 2.6- 10 15 cm). This result implies that the 
probability models do not differ significantly if we change the mass of microlenses between 17 
and 0.1M Q Thus, microlensing statistics are insensitive to changes of mass in the expected 
range of stellar masses. 



4. Discussion and Conclusions 

In the previous sections we have extended to the extragalactic domain the local (Milky 
Way, LMC, and M31) use of microlensing to probe the properties of the galaxy halos. Al- 
though our primary aim was to explore the practical application of the proposed method, 
we found that the data available in the literature can be consistently interpreted only 
under the hypothesis of a very low mass fraction in microlenses; at 90% confidence: 
a(\ogL max ) = 0.05^03 (maximum likelihood estimate) for a quasar continuum source of 
intrinsic size r s = 2.6 • 10 15 cm. This result arises directly from the shape of the histogram of 
microlensing magnifications, with a maximum of events close to no magnification and stands 
for a wide variety of microlensing models statistically representative of the considered image 
pairs. There is a dependence of the estimate of a on the source size but high values of the 
mass fraction (a > 0.2) are possible only for unexpectedly large source sizes (r s > 4-10 16 cm) . 



The low mass fraction is in good agreement with the results of EROS (ITisserand et al.ll2007l ) 
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for the Milky Way and with the limit established by MEGA fide Jong et al.l 120061 ) for M31. 
The agreement is also good with the few mic rolensing-based estimates available for individual 
object s. In RXJ 1131-1231 15ai et all (bo09h found a ~ 0.1. In PG 1115+080 [Morgan et al. 
f jgoogh obtained values in the range a = 0.08 to 0.15. For the same system iPooley et al.l 
(120091 ) found a ~ 0.1 for a source of size r, = 1.3 • 10 



cm. 



On the other hand, our estimate of the fraction of mass in microlen ses, ajlog L mn _ x ) 



0.05±n'o3, approximates the expectations for the fraction of visible matter. iJiang fc Kochanek 



( 120071 ). for instance, comparing the mass inside the Einstein ring in 22 gravitational lens 
galaxies with the mass needed to produce the observed velocity dispersion, inferred average 
stellar mass fractions of 0.026 ± 0.006 (neglecting a diabatic compression) an d 0.056 ± 0.011 



(including adiabatic compression). As discussed in IJiang fc Kochanek! (120071 ) these values 
are also in agreement with oth er estimates of the stellar mass f raction that relied on stel- 
lar popul ation models: ~ 0.08 dLintott et al.l l2006h. 0.065±o:oos faoekstra et~al1 l2005h . and 



0.031 

o!oi ( jMandelbaum et al.ll2006l ). Thus, we can conclude that microlensing is probably 



caused by stars in the lens galaxy, and that there is no statistical evidence for MACHOS in 
the halos of the 20 galaxies of the sample we considered. 

How robust are these results? There are several sources of uncertainty to consider. 
Firstly, we neglected in Equation H] the term arising from source variability, Am™". From 
a group of 17 gravitational lenses with photometric monitoring available in the literature 
we estimate an average gradient of variability of 0.1 magyear -1 . Taking into account that 
the average delay between images is about 3 months (a conservative estimate; note that the 
group of lens systems used includes many doubles, some of them with very large time delays) 
we can expect an amplitude related to intrinsic source variability of Am™" ~ 0.03, which, 
according to the histogram of magnifications (Figured]), is not significant. Moreover, if we 
assume that the probability of Am™ is normally distributed, the global effect of source 
variability is to broaden the histogram of microlensing magnifications, diminishing the peak 
and enhancing the wing. In other words, source variability leads to an overestimate of a. 
Thus, the mass fraction should be even lower if significant source variability were hidden 
in the data. In the same way, other sources of error in the data, such as the difficulty 
in separating line and continuum or in removing from the narrow emission lines the high 
ionization broad emission lines that could be partially affected by microlensing, probably 
tend to induce additional magnitude differences, Am, between the images and, hence, to 
an overestimate of a. On the contrary, cross-contamination between the spectra of a pair 
of images masks the impact of microlensing and may affect our results. Although most 
of the bibliographic sources of microlensing measurements analyze this problem concluding 
that the spatial resolution was sufficient to extract the spectra without contamination, it is 
clear that high S/N data obtained in subarcsecond seeing conditions will help to control this 
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important issue. 

Another point to address is the treatment of some of the quads, where only a subset 
of the images are used. Are we systematically excluding faint images that might be highly 
demagnified by microlensing? Let us examine the four incomplete quads in our sample. The 
fold lens SDSS J1004+4112 has two close images A and B. A is probably a saddle-point 
image and shows the most anomalous flux (Ota et al. 2006). In contrast, the optical/X-ray 
flux ratios of C and D are almost the same. Thus, there is no reason to suppose that the 
image without a useful spectrum (D) has higher microlensing probability than the others. 
PG 1115+080 is another fold quad. A\ and A2 are the two images closest to the critical curve 
and have a (moderately) anomalous flux ratio and optical variability (Pooley et al. 2007). 
The two images without available spectra (C and D) show only a small optical variability and 
are not particularly prone to microlensing. In RXS J1131-1231 the most anomalous flux ratio 
is B/C and A is a saddle-point image (Sluse et al. 2006). Image D (the one with no available 
spectrum) also has an anomalous flux but is not more susceptible to microlensing than the 
other images. Thus, in three of the four incomplete quads there is no reason to suppose that 
we are biasing the sample towards image pairs with lower microlensing probability. The case 
of SDSS 0924+0219 is more problematic. There are two sets of data for this object, one by 
Eigenbrod et al. (2006) based on observations of the low ionization lines Mgll and CIII], 
which, after two epochs of observation, reveals no difference between the line and continuum 
flux ratios of components A and B. The other set of data (Keeton et al. 2006) is based on 
Lja observations (a high ionization emission line suposed to come from a smaller region than 
the low ionization emission lines) and microlensing is detected not only in the continuum but 
also in the emission lines. This implies that the baseline for no microlensing magnification 
cannot be defined and, consequently, we could not consider Keeton et al. (2006) results. 
Anyway, we have repeated (as a test) the entire maximum likelihood estimate procedure to 
derive a but now using for SDSS 0924+0219 the microlensing measurements by Keeton et 
al. (2006). The results are almost identical: ol — 0.05^Q3. 

The size of the sample also limits the statistical interpretation. An improvement in 
the S/N of the histogram of microlensing magnifications is very important to ascertain the 
statistical significance of the low frequency of events at large magnification (only two events 
of high magnification are detected), which can impose severe constraints on the microlensing 
models. Another reason to increase the size of the sample is the possibility to define sub- 
samples at different galactocentric distances where different ratios of visible to dark matter 
are expected. In the same way it would be possible to define subsamples according to the 
type of lens galaxy or other interesting properties of lens systems. 

In any case, the impact of the main result of our study —absence of MACHOS in the 
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10 to O.1M mass range in the halos of lens galaxies— and its future prospects, points to 
the need to improve the statistical analysis in two ways: increasing the number, quality and 
homogeneity of the microlensing magnification measurements from new observations, and 
reducing the uncertainties in the macro-lens models. 
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Fig. 1. — Histogram of microlensing magnifications for the sample of image pairs in Table [3] 
(bin=0.25). 




Fig. 2. — Example of magnification maps for the case k = 7 = 0.45. From top to bottom 
and from left to right, maps correspond to a = 0.01, 0.05, 0.10, 0.15, 0.20, 0.25, 0.30, 0.50, 
1.00. 
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Fig. 3. — Relative frequency of microlensing magnifications, / Q/t)K , 7 , for pixel size (solid line) 
and r s = 2.6 • 10 15 cm (dashed line) sources (see text). Examples for three different values of 
k in the case k — j, see text. 
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HE 0047 HE 0435 HE 0435 HE 0435 HE 0512 SDSS 0006 SBS 0909 SDSS J09Z4 FBQ 0951 QSO 0957 

(B-A) (B-A) (C-A) (D-A) (B-A) (B-A) (B-A) (B-A) (B-A) (B-A) 




-2 -1 0+1 -10 +1 -1 0+1 -10 +1 -1 0+1 -10 +1 -1 0+1 -10 +1 -1 0+1 -10 +1 + 2 
Ml = +7-14 fil= + 7.17 ^1 = + 7.17 ,ul = +7.17 ,ul = -5.00 /tl = +3.33 jiil = -1.99 l^l = + Z1.05 mL = +2.59 Ail = + 6.16 
/i2 = -5.21 /i2 = -B.50 ,u2 = + 7.17 ,u2 = -4.63 /i2 = +5.94 /i2 = -1.43 /i2=+3.63 /j2 = +6.65 /i2=-0.64 /i2=-0.29 



Fig. 4. — Probability models, / aK i,oK 2 ,Ki,«;2,7i,72('^ m = m i~ m 2), corresponding to each image 
pair in the sample for different values of the fraction of mass in compact objects, a (see text). 
Iii and ji 2 are the magnifications of the images considered in each pair. The vertical dashed 
line corresponds to the microlensing measurement value. 
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SDSS J1001 SDSS J1004 SDSS J 1004 QS0 1017 HE 1104 PG1115 RXSJ1131 RXSJ1131 SDSS J1206 SDSS J1353 

(B-A) (B-A) (C-A) (B-A) (B-A) (A2-A1) (A-B) (C-B) (A-B) (A-B) 




-2 -1 0+1 -10 +1 -1 0+1 -10 +1 -1 0+1 -10 +1 -1 0+1 -10 +1 -1 0+1 -10 +1 +2 
Ml = +5.60 All = -12.87 pil = -12.87 ,ttl = +2.55 = -2.17 /il = + 15.04 ,til = -28.99 l^l = +15A3 ,uL = +4.43 ,ul = + 2.94 
/i2 = -1.25 /i2 = +25.00 ,u2 = + 3.63 ,u2 = -0.86 /i2 = +4.17 /Lt2=-16.03 /i2= + 15.43 /Li2=+11.57 /i2=-3.37 /i2=-0.88 



Fig. 5. — Probability models, / aK i,oK 2 ,Ki,«;2,7i,72('^ m = m i~ m 2), corresponding to each image 
pair in the sample for different values of the fraction of mass in compact objects, a (see text). 
Iii and ji 2 are the magnifications of the images considered in each pair. The vertical dashed 
line corresponds to the microlensing measurement value. 
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HE 1413 HE 1413 HE 1413 B J1422 B J1422 B J1422 SBS 1520 WFI J2033 WFI J2033 

(B-A) (C-A) (D-A) (A-B) (C-B) (D-B) (B-A) (B-C) (A2-A1) 




-2 -1 0+1 -10 +1 -1 0+1 -10 +1 -1 0+1 -10 +1 -1 0+1 -10 +1 -1 +1 + 2 
Ml = -5.30 fil=-5.30 ^l = -5.30 ,ttl = +9.66 /il = -15.75 /il = -0.34 ,ltl = +2.83 fll=+3.11 mL = +8.45 
/i2 = +4.7B /i2 = + 5.91 ,u2 = -3.34 ,u2 = -15.8 /i2 = +5.58 /i2 = -15.8 /i2=-0.86 /;2 = -2.63 /i2=-4.55 



Fig. 6. — Probability models, / aK i,oK 2 ,Ki,«;2,7i,72('^ m = m i~ m 2), corresponding to each image 
pair in the sample for different values of the fraction of mass in compact objects, a (see text). 
Iii and ji 2 are the magnifications of the images considered in each pair. The vertical dashed 
line corresponds to the microlensing measurement value. 
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Fig. 7. — Likelihood function vs. fraction of mass in compact objects. Dashed line (cir- 
cles): likelihood function supposing that the microlensing measurements are unaffected by 
errors. Continuous line (filled circles): likelihood function supposing that the microlensing 
measurements are affected by 0.2 mag errors. See text for details. 
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Fig. 8. — Likelihood functions corresponding to sources of deprojected size parameter, r So : 
0.8 • 10 15 cm (squares; dashed line), 3.9 • 10 15 cm (circles; continuous line), 12.4 • 10 15 cm 
(triangles; dotted line), and 40.4 • 10 15 cm (pentagons; dot-dashed line). 
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Fig. 9. — Maximum likelihood estimates of a versus deprojected source size parameter, r 
Error bars correspond to 90% confidence intervals. 



-26- 



Table 1. Micro lensing Data 



vjDject 


Image 


*** 


a ^cont .lines 


Lya 


Oil V \f \Jx V J 


V 




l\/TrrTT 

JV1 gl 1 






pair 




cont /lines * * 


cont /line* 


cont /line* 


cont /line* 


cont/linc* 


cont /line* 


cont /line* 


HE 0047-1756 1 


B - 


A 


-0.19 


1.17/1.36 














HE 0435-1223 2 


B - 


A 


-0.24 


















C- 


A 


-0.30 


















D - 


A 


0.09 
















HE 0512-3329 3 


B - 


A 


-0.40 ±0.16 
















SDSS 0806+2006 4 


B - 


A 


-0.47 ±0.20 










0.06/0.33 


0.06/0.72 




SBS 0909+532 5 


B - 


A 


-0.60 ±0.15 
















SDSS J0924+0219 e 


B - 


A 


0.00 




_ 


_ 


_ 






_ 


FBQ 0951+2635 7 


B - 


A 


-0.69 ± 0.35 










0.08/1.12 


0.13/0.46 




QSO 0957+561 8 


B - 


A 


-0.30 


-0.30/0.0 














SDSS J1001+5027 9 


B - 


A 


0.23 ±0.04 








0.63/0.35 


0.38/0.19 






SDSS J1004+4112 10 


B - 


A 


0.00 


0.50/0.50 
















C - 


A 


0.45 


0.64/0.19 














QSO 1017-207 11 


B - 


A 


-0.26 ± 0.11 




-2.21/-2.08 


-2.24/-2.06 


-2.24/-1.41 


-2.15/-1.76 






HE 1104-1805 12 


B - 


A 


0.60 ±0.03 


— 








1.75/1.12 


1.68/1.12 


— 


PG 1115+080 13 


A2 - 


- Al 


-0.65 


-0.65/0.0 














RXS J1131-1231 14 


A - 


B 


1.39 














0.65/-0.74 




C- 


B 


1.58 














1.27/-0.31 


SDSS J1206+4332 15 


A — 


B 


-0.56 ±0.21 










0.32/1.08 


0.54/0.89 




SDSS J1353+1138 16 


A — 


B 


0.00 
















HE 1413+117 17 


B - 


A 


0.00 ±0.04 






0.23/0.19 


0.20/0.23 










C- 


A 


-0.25 ±0.10 






-0.03/0.27 


-0.07/0.27 










D - 


A 


-0.75 ±0.08 






0.2/-1.07 


0.22/-0.85 








B J1422+231 18 


A — 


B 


0.16 




0.27/0.11 














C - 


B 


0.02 




0.75/0.77 














D - 


B 


-0.08 




3.92/4.00 












SBS 1520+530 19 


B - 


A 


-0.39 ±0.07 








-0.04/0.27 


0.08/0.54 






WFI J2033-4723 20 


B - 


C 


-0.50 








-0.09/0.41 










A2 - 


- Al 


0.00 








0.32/0.32 









* Magnitude differences between images in the continuum and in the line emission, respectively (when an individual value for one or more 
lines is available) 

** Magnitude differences between images in the continuum and in the line emission, respectively (when a global value for an spectral 
region including several lines is given) 

""Average microlensing magnification, < Am >=< Am™" 1 — Am Ime > 

IWisotzkT et al. | J2004h (Line flux ratio given by the authors. Continuum flux ratio estimated from Figure 3) 

^Wisotzki et al] l l2003t) (Microlensing magnifications taken from Table 3) 

'Wucknitz et al.l ll2003T ) (Microlensing magnification estimated from Figure 3) 

ilnada et al.l ll2006T ) (Flux ratios computed from electronically digitized spectra) 

;: Mcdiavill a et al] ll2005h (Microlensing magnification estimated from Figure 7) 

teigenbrod et al. | <2006ft (See text) 

' Schcchtcr et al. ( 1998) (Flux ratios computed from electronically digitized spectra) 
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^Goicoechea et aLl i2005l) (Line flux ratio given by the authors. Continuum ratio estimated from Figure 1) 
^Oguri et al,l d2005T) (Flux ratios computed from electronically digitized spectra) 
lc lG6mez- Alvarez et al.l d2006T) (Flux ratios estimated from Figures 3 and 4) 
n Surdcj et al.l l ll997l) (Flux ratios computed from electronically available spectra) 

1: ^Wisotzki et all lll993h (Line flux ratio given by the authors. Continuum flux ratios estimated from Figure 3) 
ls |Popovic &: Chartasl j2005h (Line flux ratio given by the authors. Continuum flux ratio estimated from Figure 9) 
HSluse et al] [120071) (Flux ratios taken from Table 5) 

1 ^Oeuri et alj J2005T) (Flux ratios computed from electronically digitized spectra) 
le llnada et all J2006h (See text) 

1 1Popovic &: Chartasl (|2005h (Flux ratios computed from electronically available spectra) 
1 ^Impev et al. (1996) (Flux ratios taken from Table 3) 

1 ^Chavushvan et alj lll997h (Flux ratios computed from electronically digitized spectra) 
2( lMorgan et al.l II2004TI (Flux ratios estimated from Figure 9) 
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Table 2. Comparison between Emission Line and mid-IR Flux Ratios 



Object 


Image Pair 


Am linea 


/\ ni mid— I R 


SDSS J1004+4112 1 


B - 


A 


0.50 


0.30 




C - 


A 


0.19 


0.50 


HE 1104-1805 2 


B - 


A 


1.12 


1.13 ±0.06 


PG 1115+080 3 


A2- 


-Al 


0.0 


0.08 ±0.06 


HE 1413+117 4 


B - 


A 


0.21 ± 0.02 


0.19 ±0.07 




C - 


A 


0.27 ±0.00 


0.36 ±0.07 




D - 


A 


0.96 ± 0.11 


0.99 ±0.06 


B J1422+231 3 


A — 


B 


0.11 


0.18 ±0.05 




C - 


B 


0.77 


0.61 ±0.06 



1 Mid-IR data from lRoss et alj l l2009h 
2 Mid-IR data from IPoindexter et al] d2007h 
3 Mid-IR data from lChiba et al] l l2005h 
4 Mid-IR data from lMacLeod et al] (2009) 
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Table 3. Lens Models 



Object Ratio Value R innf , r (kpc)* -Router(kpc)* k\ k 2 71 72 



HE 0047-1756 




-A 


-0.19 


TB 


= 


3.6 


r A 




4.9 


KA 


= 


0.43 


k b 


— 


0.61 


7A 


= 


0.48 


7s 


= 


0.65 


HE 0435-1223 


Amj- 


-A 


-0.24 


tb 


= 


7.6 


T A 




8.2 


KA 


= 


0.46 


k b 


— 


0.52 


7A 


= 


0.39 


7s 


= 


0.59 




Amc- 


■ A 


-0.30 


r c 


= 


8.2 


TA 




8.2 


KA 


= 


0.46 


KC 


— 


0.46 


7A 


= 


0.39 


7C 


= 


0.39 




Amo- 


-A 


0.09 


TD 


— 


7.0 


TA 




8.2 


KA 


— 


0.46 


KD 


— 


0.56 


7A 


— 


0.39 


7d 


— 


0.64 


HE 0512-3329 


Amg- 


■A 


-0.40 


r A 




2.2 


TB 




3.3 


KA 


= 


0.59 


K B 


= 


0.41 


7A 


= 


0.55 


7b 


= 


0.37 


SDSS 0806+2006 


Am B - 


■A 


-0.47 


r A 




7.2 


TB 




3.0 


KA 




0.34 


KB 




0.82 


7A 




0.29 


7b 




0.77 


SBS 0909+532 


Am B _ 


-A 


-0.60 


r A 




3.4 


ra 




6.3 


K A 




1.00 


K B 




0.55 


7A 




1.39 


7b 




1.19 


SDSS J0924+0219 


Am B - 


■A 


0.00 


r A 




5.2 


ra 




5.8 


K A 




0.50 


K B 




0.45 


7A 




0.45 


7b 




0.39 


FBQ 0951+2635 


Atrifl- 


■A 


-0.69 


r A 




3.6 


TB 




0.9 


KA 




0.28 


k b 




1.07 


7A 




0.15 


7b 




1.02 


QSO 0957+561 


Am B _ 


■A 


-0.30 


rs 




5.7 


rA 




28.6 


K A 




0.20 


k b 




1.03 


7A 




0.15 


7b 




0.91 


SDSS J1001+5027 


ArriB- 


■A 


0.23 


r A 




10.6 


TB 




5.0 


KA 




0.35 


K B 




0.74 


7A 




0.28 


7b 




0.72 


SDSS J1004+4112 


Am B - 


■A 


0.00 














KA 




0.48 


kb 




0.48 


7A 




0.59 


7b 




0.48 




Amc- 


■A 


0.45 














K A 




0.48 


KC 




0.38 


7A 




0.59 


7C 




0.33 


QSO 1017-207 


ArriB- 


-A 


-0.26 


r A 




5.4 


TB 




1.5 


K A 




0.35 


K B 




1.23 


7A 




0.45 


7B 




1.32 


HE 1104-1805 


Atnj_ 


■A 


0.60 


r A 




8.6 


ra 




16.6 


K A 




0.64 


k b 




0.33 


7A 




0.52 


7b 




0.21 


PG 1115+080 


Am A2 


-Al 


-0.65 


r A i 




= 5.8 


rA2 = 


= 5.9 


Ml 




: 0.49 


KA2 




: 0.51 


7A1 




: 0.44 


7A2 




: 0.55 


RXS J1131-1231 


ArriA- 


B 


1.39 


r A 




10.1 


rB 




10.2 


K A 




0.44 


K B 




0.43 


7A 




0.59 


7B 




0.51 




Am C - 


B 


1.58 


r A 




10.1 


rc 




10.7 


k b 




0.43 


*C 




0.42 


7s 




0.51 


7C 




0.50 


SDSS J1206+4332 


AmA- 


B 


-0.56 


TB 




10.02 


rA 




15.0 


KA 




0.43 


k b 




0.63 


7A 




0.41 


7B 




0.72 


SDSS J1353+1138 


ArriA- 


B 


0.00 


TB 




1.6 


rA 




5.2 


K A 




0.30 


k b 




0.96 


7A 




0.22 


7B 




0.89 


HE 1413+117** 


Amg_ 


-A 


0.00 














K A 




0.53 


k b 




0.43 


7A 




0.64 


7b 




0.34 




Amc- 


-A 


-0.25 














K A 




0.53 


KC 




0.46 


7A 




0.64 


7C 




0.35 




Am,£j_ 


-A 


-0.75 














K A 




0.53 


K D 




0.58 


7A 




0.64 


Id 




0.69 


B J1422+231 


Arrij\- 


B 


0.16 


I'B 




5.2 


rA 




5.3 


K A 




0.38 


k b 




0.39 


7A 




0.53 


7b 




0.66 




Amc- 


B 


0.02 


TB 




5.2 


rc 




5.7 


K B 




0.39 


KC 




0.36 


7s 




0.66 


7C 




0.48 




Amj)_ 


-B 


-0.08 


VB 




5.2 


ro 




1.3 


K D 




1.54 


k b 




0.39 


Id 




1.81 


7b 




0.66 


SBS 1520+530 


ArriB- 


-A 


-0.39 


r A 




9.5 


rB 




3.0 


K A 




0.29 


kb 




0.90 


7A 




0.15 


7b 




0.85 


WFI J2033-4723 


Atnj_ 


■C 


-0.50 


rc 




6.8 


rs 




11.2 


K B 




0.38 


kc 




0.61 


7s 




0.25 


1C 




0.73 




Am A2 


-Al 


0.00 


r A 2 = 


= 8.2 


rAi 




= 9.3 


KA1 




: 0.48 


KA2 




: 0.55 


7A1 




: 0.39 


7A2 




: 0.65 



* Computed using a concordance cosmology 
** Lens redshift unknown 
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Table 4. Frequency Distribution of Image Pairs, n fcljfc2 



k 2 = 0.35 k 2 = 0.45 k 2 = 0.55 fc 2 = 0.75 k 2 = 0.85 k 2 = 1.05 k 2 = 1.15 k 2 = 1.25 fc 2 = 1.55 k 2 = 1.75 



fci = 0.35 2 1 3 2 3 1 1 

fci = 0.45 4 9 1 

fci = 0.55 1 1 



